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NUMERICAL SIMULATION OF STEADY AND UNSTEADY VISCOUS FLOW IN 
TURBOMACHINERY USING PRESSURE BASED ALGORITHM 


B. LAKSHMINARAYANA, Y. HO and A. BASSON 




The Pennsylvania State University 
Department of Aerospace Engineering 
University Park, PA 16802 



The objective of this research is to simulate steady and unsteady viscous flows, including 
rotor/stator interaction and tip clearance effects in turbomachinery. 

The numerical formulation for steady flow developed here includes an efficient grid 
generation scheme, particularly suited to computational grids for the analysis of turbulent 
turbomachinery flows and tip clearance flows, and a semi-implicit, pressure-based computational 
fluid dynamics scheme that directly includes artificial dissipation, and is applicable to both viscous 
and inviscid flows. The values of these artificial dissipation is optimized to achieve accuracy and 
convergency in the solution. The numerical model is used to investigate the structure of tip clearance 
flows in a turbine nozzle. The structure of leakage flow is captured accurately, including blade-to- 
blade variation of all three velocity components, pitch and yaw angles, losses and blade static 
pressures in the tip clearance region. The simulation also includes evaluation of such quantities of 
leakage mass flow, vortex strength, losses, dominant leakage flow regions and the spanwise extent 
affected by the leakage flow. It is demonstrated, through optimization of grid size and artificial 
dissipation, that the tip clearance flow field can be captured accurately. 

The above numerical formulation was modified to incorporate time accurate solutions. An 
inner loop iteration scheme is used at each time step to account for the non-linear effects. The 
computation of unsteady flow through a flat plate cascade subjected to a transverse gust revels that 
the choice of grid spacing and the amount of artificial dissipation is critical for accurate prediction 
of unsteady phenomena. The rotor-stator interaction problem is simulated by starting the 
computation upstream of the stator, and the upstream rotor wake is specified from the experimental 
data. The results show that the stator potential effects have appreciable influence on the upstream 
rotor wake. The predicted unsteady wake profiles are compared with the available experimental data 
and the agreement is good, the numerical results are interpreted to draw conclusions on the unsteady 
wake transport mechanism in the blade passage. 
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OBJECTIVE 
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NUMERICAL METHOD 
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• Validate for inviscid flow computation 
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NUMERICAL METHOD (CONTD.) 
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VALIDATION/TEST CASES 
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e compared to the flat plate cascade 




VALIDATION/TEST CASES (CONTD.) 
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Tip clearance flow through a turbine cascade 




VALIDATION/TEST CASES (CONTD.) 

Unsteady response function of a flat plate cascade and a C4 cascade 
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Fraction of Axial Chord 
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UNSTEADY VISCOUS OVER AN AIRFOIL 

The unsteadiness is generated by two flapping foils oscillating about 
chord upstream the tested foil 
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UNSTEADY VISCOUS OVER AN AIRFOIL 
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X/C=0.389 X/C=0.61 1 X/C=0.784 

PRESSURE SIDE 








AMPLITUDE OF VELOCITY FOR HARMONIC n= 
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TIP CLEARANCE MODELLING 
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3D STEADY FLOW IN THE TIP CLEARANCE REGION 

FLOW OF A TURBINE 
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3D STEADY FLOW IN THE TIP CLEARANCE REGION 

FLOW OF A TURBINE 
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good tip gap resolution 



EMBEDDED H-GRID TOPOLOGY 


Grid Below Tip Gap 
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COMPARISON OF MEASURED AND PREDICTED 
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COMPARISON OF MEASURED AND PREDICTED 
PITCHWISE FLOW ANGLES 




Pitchwise Flow Angle (deg) 
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Numerical 
Bindon (1990) 
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Numerical 
Bindon (1990) 




COMPARISON OF MEASURED AND COMPUTED SPANWISE FLO 
ANGLES (POSITIVE ANGLE WHEN FLOW DIRECTED TOWARD 

ENDWALL) 
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Numerical 
Bindon (1990) 



CONTOURS OF LOSS COEFFICIENT (cj 
PREDICTED AND MEASURED 
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Pitchwise Pjtchwise 



DEVELOPMENT OF MASS AVERAGED LOSS 
COEFFICIENT (c t ) 
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Axial Chord 
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OF POOR QUALITY 







VELOCITY VECTORS AT CONSTANT SPANWISE 
POSITIONS PROJECTED ON TO A BLADE-TO-BLADE 

(XY) PLANE 
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VELOCITY VECTORS AT CONSTANT SPANWISE 
POSITIONS PROJECTED ON TO A BLADE-TO-BLADE 
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VELOCITY VECTORS AT CONSTANT SPANWISE 
POSITIONS PROJECTED ON TO A BLADE-TO-BLADE 

(XY) PLANE 
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Midspan Exit Velocity 







COMPUTED SECONDARY YORTICITY VECTORS 

[ ( VX ^~^mld s pa»)) ] 






COMPUTED SECONDARY VORTICITY VECTORS 
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20000. per sec 




CONCLUSIONS 

• Geometric series distribution scheme 

Good control over the boundary points 

• Algebraic-elliptic grid generation scheme 

Good control over clustering 
Good control over orthogonality 
Enhanced stability of elliptic generation 

• Embedded H-grid 

Single block discretization for tip clearance cases 
No Modification of blade tip shape required 
Retains H-grid connectivity pattern 

• Effect of artificial dissipation 

Numerical accuracy vs. convergence rate 
Minimum artificial dissipation should be used 
• Modelling of tip clearance flows 

Major physical phenomena captured 

Location of leakage vortex, pltchwise and spanwise angles, losses, 
static pressures predicted accurately 
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• 2D unsteady flow computation 
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velocity profiles and blade pressure agree with the measurements, 
including the magnitude and phase angle. 




